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Abstract 

The topology of the surface electronic states is detected with photoemission. We explain the 
photoemission from the topological surface state . This is done by identifying the effective coupling 
between surface electrons-photons and vacuum electrons. The effective electron photon coupling is 
given by er^ where r is the dimensionless tunneling amplitude of the zero mode surface states to 
tunnel into the vacuum. We compute the polarization and intensity of the emitted photoelectrons. 
We introduce a model which takes in account the Dirac Hamiltonian for the surface electron to 
photons coupling and the tunneling of the zero mode into the vacuum. Within the Green’s function 
formalism we obtain exact results for the emitted Photoelectrons to second order in the laser field. 
The number of the emitted photoelectrons is sensitive to the laser coherent state intensity, the 
polarization is sensitive to the surface topology of the electronic states and the incoming photon 
polarization. The calcnlation is performed for the helical , Zeeman and warping case allowing to 
study spin textures. 
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I. Introduction 


Photoemission, photoconductivity, optical conductivity and scanning tunneling mi¬ 
croscopy are sensitive to the nature of surface states. Photoemission is studied using a 
high power laser-based light source. It has been shown that the spin polarization of the 
photoelectrons emitted from the surface of 612803 topological insulator (TJ) P-|l] can be 
manipulated through the laser light polarization UM finds that the photoelectron polar¬ 
ization is completely different from the initial state and is controlledr by the photon polar¬ 
ization. Few explanation based on phenomenological models have been proposed [3^3- In 
a recent photoemission experiment ra the authors have demonstrated ± 100 % reversal 
of a single component of the measured spin polarization vector upon the rotation of light 
polarization as well as full three-dimensional manipulation by varying experimental configu¬ 
ration and photon energy. This experiment shows that the photoelectrons spin polarization 
is achievable in systems with a layer-dependent, entangled spin-orbital texture. There are 
also other studies of spin-polarized photoelectrons spectroscopy of TI m including orbital- 
selective spin textures [ 12 ] and reports that interactions [13] might affect the photoemission 
spectrum. Regarding the helicity-dependent photocurrent, due to the spin selection rules 
one hnds that circularly polarized light excites the surface states and an electric DC current 
is observed [H] and was investigated by [15] . The authors compute the induced DC current 
using the TI surface model for 612803 which also includes the warping nonlinearity and the 
presence of the Zeeman magnetic term [T3] . 

In spite of this success the main question of how to explains photoemmission for Topo¬ 
logical Insulators TI remains open. It is not clear what is the electron-photon coupling 
which describes the coupling between the electrons in the vacuum and in the solid. On 
the surface of the TI the coupling is given by the Dirac form a ■ A and in the vacuum 
hj A ■ p, as a result neither Hamiltonians can describes the transition between the sur¬ 
face and vacuum states. Physically one describes photoemmision as a process where a 
photon is absorbed and an electron is excited from the surface to the vacuum, clearly no 
such matrix element exists for the surface electrons (for the bulk electrons the situation 
is different since the electron photon coupling is given hy A ■ p and the matrix element 
{vacuum — electronA ■ pbulk — electron) 7 ^ 0). 8 uch a process require the knowledge of 
the electron photon vertex, surface{electron) — photon — vacuum{electron). We solve this 
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problem using the zero mode surface TI state. The surface state are localized at z = 0 with 
the amplitude [12] (^ < 0 represents the solid and z > 0 describes the vacuum, the 
eigenstates on the surface have a small amplitude to tunnel into the vacuum.) z = 0 rep¬ 
resents the surface and z > d describes the free electrons separated by the surface-vacuum 
binding energy Vq. 

The effective coupling is given by e • (e is the charge and r is the overlap between the 
two types of wave functions). 

We hnd that the polarization measured by the detector depends on the product : the 
projection of the electron spin polarization on the direction of the detector, the scalar product 
between the photon vertices (which result from the spinor form of the electron operator) 
and the transverse polarization of the incoming photon) and the intensity which measures 
the number of the photoelectrons emitted. 

The plan of this paper is as follows. In Sec. II we present the model for surface in 
the presence of the photon held and tunneling amplitude into the vacuum. In Sec. Ill we 
introduce the Green’s function and compute the number of the photoelectrons emitted for 
the helical, Zeeman and warping case. 

Section IV deals with the study of the polarization for the helical, Zeeman and warping 
case as function of the polarization and intensity of the incoming photons. Section V is 
devoted to discussions and presents our main conclusions. 

II- Model for the photoemission from a TI surface 

The photoemission for the TI involves a four component spinor for the bulk, a two com¬ 
ponent spinor for the surface and a wave function for the vacuum. The surface electrons 
at z = 0 have an amplitude I to tunnel into the vacuum (z > 0). The electrons 

^ y/NiK,Ko) ^ ^ 

detected by the detector have a mean free path of L. To simplify the problem we expand 
in plane waves the vacuum electrons using a box of length L. We have for the vacuum 
electron, K{K,z > 0) = Efc, a{K,kz). The electrons on the surface are described by 
the spinor ^(rT) with the eigenspinor \u{K)), ^(hl) = / u{K) is the 

two component spinor Ua=^{K) = (■!■ \u{K)) and Ua=i{K) = (j, \u{K)). For z > 0 we have 
'^{r±,z) ^ z = 0) and for 2; < 0 we have T(r_L,^) -^£==^[f±, z = 0). 

The surface electrons overlap the with the vacuum electrons in region closed to z=0. 
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We restrict the overlap to the region —do < z < d where do = -^ and d = ^ and find the 
tunneling matrix element t{kz)- (;^ inverted gap of the TI and ^ is gap of the vacuum and 
d, do are a few lattice constants). 



dzC^K, z)u*^{K)ba{K, z) + h.c. 


= C\K,z = 

= C\K,z = 


0)<(i?) 


rd=h 


dz 


^—K,{z)z pikzZ 


l-do=- — 

KQ 


vL k. 


h.c. 


0)uUK)J2t{kz)b^{K, kz) + h.c. 


( 1 ) 


where the tunneling matrix is t{kz) is given by , 


t{kz) = / " , dz—. 


^—k{z)z ^ikzZ 


( 2 ) 


In{k, ko) 

t{kz) is controlled by k(; 2 ), n{z > 0) = k, k{z < 0) = kq and the normalization factor 


^ N{k,ko) 

Using Eq.(l) we introduce the following model: 


— ffisurf.) _|_ jj{vac.) _|_ jj(surf.-vac.) _|_ 

r d'^K ^ ^ ^ 

= J j^^\K,z = 0)[hv{a^K2-(J^Ki) + a^A{K,kc)]^{K,z = 0) 

r d'^K ^ 

^ j C't(A',2 = 0) y U:{K)Y.t(kOK(KX) + h.c] 

' 5'i- kz 

//<“'■>(«) = (-ev) I ^ I ^C\K + Q, 2 = 0; t)C(K, z = 0; t)WiK + Q, A) ■ A(Q, () 


( 3 ) 


ffisurf.) (describes the TI surface electrons. We will consider three different cases: 

A- The helical state A(K,kc) = 0 is zero, the eigenvalues are given by e = Cgurf. = 
Co = v\K\, e = eo — /i = v\K\ — jj, and the spinors for \K\ ^ 0 are }u{K)) = -^\K) ® 
[1, —. The electromagnetic vertex functions are given by the Pauli 
matrix elements, Wi{K + Q, K) = (u{{K + Q)\{—a‘^)\u{K)) = cos[^(x(i? + (5)+ x(i?))] and 
W 2 {K + g, A) = {u{K + g)|aXA)) = - sin[|(x(A + g) + x(A))]. 
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B- For the Zeeman gap A 7 ^ 0 the eigenvalue are given by e = eo + and the 


spinors are, |m(A, A)) = \k) 


cos 


V 2 


COS 


PiK) 


^y{nvK)^2+A^2 ' 


V 2 

The vertex functions are for this case, Wi{K,K) = sin[/3(iF)] sin[x(iF)] , W 2 {K,K) = 
-sm[/3{K)]cos[x{K)]. 


C- For the nonlinear warping the eigenvalue are. 


e = 


Co 


^ + (t)^cos2[3x] 

(cc = vkc is the warping energy.) The spinors are given by, {K, kc)) = \K) 


cos 




V 


2 


V 


The vertex function take the form, Wi{K,K) = 


sin[/3(iF,/Cc)] sin[x(iF)] , W 2 {K,K) = —sm[/3{K,kc)]cos[x{K)] with cos[P{K,kc) = 


(A)^cos[3x(fc)] 


The vacuum Hamiltonian jg controlled by the binding energy Vq and eigenvalues 

E = E — fi = + ho ~ h- The vacuum electrons operators obey the momentum 

expansions ba-{K,z > 0) = ,kz) where L is the distance from the TI surface 

to the detector. 

is the electron-photon Hamiltonian restricted to the surface aX z = 0. 
Ai{Q, Qz] t), Wi{K + K) are the photon field and vertex for the f = 1, 2 direction. 

The high intensity photon held A{Q, qz) is a coherent state IH). The direction of the incoming 
photon with respect to the surface at 2 : = 0 is given by the vector p = Q + Qz- The two 
transverse linear polarizations are given by the vectors 6^=1 (p) and es= 2 {f^ which obey 
• || = 0, e,=i(p) ■ es=2P = 0. 


f d^P 1 
A{f,t) = \-‘ 




(27r)3 y^2H(p) 

= (sin( 6 *) cos( 0 ), sin( 6 *) sin( 0 ), cos( 6 ')), 

es=i(p) = (cos( 6 ') cos( 0 ), cos( 6 ') sin( 0 ), - sin( 6 ')), es= 2 {p) = (-sin( 0 ), cos( 0 ), 0 ) 


(4) 


6, 0 are the photon polarization angles, |H) is the the lase coherent state and e is the 
dielectric constant. 


Ill- The spin detection 

The detector is in the plan x — y parallel to the TI surface and perpendicular to the 
the axis . The detector measures the spin polarization in the y direction. When the TI 
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surface rotates around the y axes the detector measures the angle (pd which coincide with 
the spinor angle x[K] = (pd- The rotation of the sample around the x axes allows to measure 
the momentum in the z direction. We have \K\ = |/c|cos[6*d] and = |fc|sin[6*d] where 
k = K + zkz- The eigenvalue of the vacuum electrons can be written in term of eo (the 
surface eigenvalue ) and 9d, E = E - y = +Vo- y. 

The spin density (n^) measured by the detector is given in terms of the Green’s function 
Gfi,a{Q,qz;t,t + 5t): 


(n^(^ 
= (-« 



/ dk’ -* 

kz-,t,t + 6t) 
kz-,u:),6t^0 


Gfi^a{K, Qz] t,t + 6t) = -i{g\T{bfd{K, kz] t)bl{K, kz] t + 6t))\g) 


= -t{0 ® n\T(bfd{K, kz; t)bUK, kz; t + 5t)e^ ^ 


(5) 


T stands for time order operator, \0) represents the ground state of the Hamiltonian, 
]^(surf.) _|_ jj{vac.) _j_ j^(surf.-vac.)^ \0)c means ”Connected” diagrams. 

Using Wick’s theorem [19] with respect the ground state \0) we compute the Green’s function 
Gj 3 ^a{K, kz] 0 j) to second order in the photon held. 

Gp,^{K,kz-,u) = — M,{i,j\e,<P)Wi{K)Wj{K) 


gG’^\K,kz\K-,u:) 


dui 




( 6 ) 


Ms(i, j|6*, 0) is the photon matrix polarization , Ms(i, j|6*, 0) = {es{9,(p) ■i){es{9,<p) ■ j). 

The Green’s function gG°9(^x\K, kz] oj) describes the creation of a vacuum electron with 
spin a and the destruction of a surface electron, gGx)(^j{^ kz\K]Uj) describes inverse process 
and gG^\K\K]Uj) represents the Green’s function for the surface electrons. Ds{uj) is the 
photon Green’s function dehned with respect the coherent state |G). 

g^^^G^K\K,kz]t) = -z{O\T{C{K]t)bi{K,kz]0))\O), 
gG>-\K,kz\K]t) = -i{O\T{b^iK,kz]t)C^{K]0))\O) 

g<^^’^\K\K]t) = -z{O\T{C{K]t)G\K]0))\O) 
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D.(t) = -*{SJ|T((o,(p)e-“ + a'.(p)e-“)(o,(p) + o;(p)))|Ji) 


( 7 ) 


The Green’s function /c^; t) kz\K;t) and g^'^’'^\K\K;t) are obtained from 

the Heisenberg equations of motion, 

dbai^K^ kzi t) 


ih- 


dt 


dC{K, z = 0]t) 

Lit I 

dt 


ba{K, kz] t), + if 

C{K, z = 0]t), {surf.-vac.) 


( 8 ) 


We obtain: 


g^'^’^\K\K,kz-,u) = 

_ t*{kz)U:{K) _ 

(^u — e + i6sgn[e]j (^oo — E + i6sgn[E]^) — E{kz) 

f a;) = 


f{kz)Ul{K)V[K\K,kz{u] 


t{kz)U^{K) 


— e + i6sgn[e]^ (^u — E + i6sgn[E]^ — E{kz 


= tikz)UaiK)T[K, kz\K;u] 


T[K,kz\K;uj]=T[K\K,kz;u = 


0[e] 


+ 


(l-0[6l) 


{u — E^ + iS){u — E_ — id) {u — E^ — i6){u — E_ — i 


( , 


where E± = ± }^\J{E — e)"^ + AE{kz) are the two roots of the algebraic equation, {oj — 

e){u — E) —t^{kz) = 0. Ua{K) is the the spin component a =t,i of the surface spinor u{K). 
0 [e] is the step function (one for e > 0 and zero other wise) which at hnite temperatures is 
replaced by the Fermi-Dirac occupation function. 

The Green’s function for the surface electrons is given by: 


= u — e + idsgnli] — 


t^{kz 


n-l 


u — E + iSsgn[E]^ 


( 10 ) 


^22 

where E = + Vo — /i. We replace the sum by a momentum integration and 

use a contour integral. The tunneling matrix element given t{kz) given in eq.(l) is replaced 
by constant tunneling matrix element , t{kz) ~ r/i, where r is dimensionless {t{kz) has 
dimension of energy). The replacement of the Eq.(9) by the integral / dkz introduces 
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C.D On 



the quantization length L. The restriction of the 2 : integration in Eq.(9) by \z\ < d will 
introduce in Eq.(lO) the dimensionless parameter ^ which will be approximated by KpL 
where Kp is the surface Fermi momentum. 


ho) (Eq p) T 


u — € + i5sgn[e] + (KpL) 


0 [ca tajyiax Eq, hi] 1 1 


^max Eq, ff] 


(e + p)^- 

2 fi 


( 11 ) 


We substitute into Eq.(5) the Green’s functions g^°''‘^\K,kz\K-,u), g^^’°‘\K\K,kz;uj), 
g^^’^\K\K;u) and perform the contour integral with respect the photon Green’s function, 
J duj,g('^’'^^{K\K-,uj-uj,)DsM, 


D..{u) = [(- 


l + iV, 


+ 


iV. 


in ) ( 


A, 


+ 


1 + iV, 


( 12 ) 


kl + ig u + kl + ig'' — kl + ig u + kl + ig^'i 
(In obtaining Ds{u) we have used the coherent states properties of the photon held with 
the occupation number Ng, aJa^lG >= Ns\Ll >.) The momentum integration in Eq.(l), 
/ kz] co) is replaced in terms of the polar angles y = 0^, 9^ and the surface 


energy e by / ^ / 


dx r dip. 


27r d 2 k \_ {Tiv)‘^ 


f 


de 


Je±G 


(€Q+^) deQ 


dA G 0 ^a{e + /i, kz;u) . We hnd from Eq. (5) that the 


spin density measured by the detector, 


{ny{9d, 0d|6', 0, s; e + p)) = P|(6', 0, 9d, (t)d, e, A, ec)dl{9d, (t)d, e + p) 
where P|(6', 0, 6*^, 0^, e. A, Cc) is the spin density polarization controlled by the photon held, 
dl{6d, <pd, e + /i) is the number of photoelectrons per solid angle for the surface energy e. The 
density polarization for the Zeeman gap P^(6', 0, 6 d, <j)d, G A) can be written as a product of 
the spin density in for the gap A = 0 (the helical state) or the warping Cc, Ps(6', 0, 9d, 0d, e + 
/i. A) = P^g{0, 0, 6 d, (j)d, A = 0) and the gap function, PpoL-supr.[e + /i, A], 


{ny{9d, 0d|6', 0, s; e + p. A)) = P^(6', 0, 9d, 0d, e, A)dl{9d, (j)d, e + fx) 

P^g{6, 0, 6d, 4>d, e + p. A) = P^g{6, 0, 6d, 4>d, A = 0, Cc = 0)PpoL-supr.[e + K A] 

Ppol.—supr. [e + p.. A] — 


(e + p) 


2 , 3 

2 


(e + g)‘^ + A2 

2 4 

(il(6'd, 0d, e, + Ns)n],[L\Q 


- ^ r- A TT / Kp . p(e + u)\/P — A^ 

E- + G — Umax)^-, A, Eo, G]] (- -)- 


( 


p 2 


0 ; 


p 


(E+-E_)3 
d(j)d dOdde 


K- + a- A, 1^,'■(g- + f! -^ 

( 13 ) 



For the warping case we have, 


(n^((9rf, (j)d\0, (j),s;e + /x)) = P^((9, 6 *^, 0^, e, ec)(il((9^, e + /x, 


Pg (^7 05 4^dl - ^5(^505 ^C?5 0C?5 O) Ppo/,_5^ip7’_ [6 “h /X, 6^] 

p ■ ■ ' 


pol.—supr. 


. e + p, Cc] — 


'1 + (^)'^cos2(30d) 

2 4 

rfl(0d,0,,e,ee) = (^)(iFFi:)(l + iV.)4[e]0 


A ^ A / -^F , uie + w)Ve^ — 

-E- + r2 — CUmaxie, A, Vo, h2] (p )- 


J (E+-E_)= 


( 




:)(; 




F + n - A. K,. n|' ' (P_ + f! - e)= + 

dcpd dOfide 


1 3(^)4 cos 2 [30^] 2n Ticos^[6d 


(14) 


ni?[e] is the Fermi Dirac occupation function and 0pj is the step function which is one for 
a; > 0 and zero otherwise, dl represents the number of emitted photoelectrons per solid angle 
dd,4>d and surface energy e. This number depends on the number of the incoming photons 
Nsi the ratio between the Fermi momentum Kp and the photon momentum Pph_ = ^ the 
effective charge e-r^ and chemical potential p. We have plot the function dl for the followings 
values: 


jj, = O.leD is the chemical potential,^ = 5eD is the binding energy D = 5eD is the laser 
frequency, r = 0.1 is the dimensionless tunneling parameter, KpL = 10^ is determined by 
the distance from the surface sample to the detector, n = 5 x 10^ — is the Fermi surface 
velocity = jA jg the hue structure constant and Ng = 10 ^^ is the number of photons. 

The plots will be restricted to energies e > —p which represents the surface electrons, the 
bulk contributions will be ignored. 

We hnd that 10^^ photons are needed for two photoelectrons to be emitted. Figure 1, 
2 , 3 and 4 shows the number of photoelectrons emitted for A = 0, A = O.OSeD, warping 
energy Cc = O.OSeD,^^ = | and = |. In all the four cases we have used the same values. 
The warping the angle (pd controls the number of the photoelectrons . Each time the angle 
(pd = (2n + 1)|, n = 0 , 1,2, .. the intensity is maximum and is minimum when the angles 
are (pd = (2n + 1)|. 


IV- The spin density in the y direction measured by the detector at the polar 
angle (pd-,0d as a function of the incoming photons 6^=1 2 (^, 0 ) 
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dI[£-^.e^=;r/4;A=0,Af^=10‘'^ ] 



FIG. 1: The number of photoelectrons dl{6d, (l>d, e) for A = 0 

dI[E-Ai,ed =i/4;A=0.05,Afs =lo'^ ] 



FIG. 2: The number of photoelectrons dl{6d-, (j^d, e) for A = O.OSey 

The spin polarization of the TI surface is given in terms of the spinor states Ua{K) 
with a =t, i- The momentum parallel to the surface is conserved, the chiral angle satisfies 
cos[0d] = cos[x]. The spin polarization (a^) = = cos[0(i](sin[/?[e, A, ec ])2 0[e] is recorded 

by the detector as 0, (j)d, e). sin[/3[e, A, ed represent the effect of the Zeeman gap A or 
,ec for warping. The function (sin[/3[e, A, ec])^ is given for the gap A and warpping Cc, 

Ppoi.-supr.[e +fi,A] = (sin[/?[e; A)])i = ^ 2 )" 

PpoL-supr.[i + f^,ed ^ (sin[/?[6,ee)])2 = (^ 

(15) 

Using the definition given in Eqs.(4 — 6) we obtain the relation between and 

A = Cc = 0): 

py{e,<P,<P,,e,A = e, = 0) = Sy ^ M,{i,j\e,<P)Wi{<Pd,e)W,{<Pd,e)W,{<Pd,e) 
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=nlA,ip^ =n-/3;ec=0.01 .Afj =10'^ ] 



FIG. 3; The number of photoelectrons dl{9d-, 4>d, e) for Cc = O.ObeV, (pd = ^ 


dlle-ti.Sj =i!l'i,ipi =n-/6;ec=0.01 .Afj =10'^ ] 



FIG. 4; The number of photoelectrons dl{6d-, (pd, e) for Cc = O.OSeG, (pd = ^ 


( 16 ) 

The spin density for the linear for the p polarization polarization corresponds to es=i{(p = 
0, 6'),and the s polarization is corresponds to es= 2 {(p = 0). We hnd: 

{9, (p = 0, (pd, A = ec = 0) = cos^ [6^] cos^ [cpd] 

P^=2(0 = O,0<i, A = Cc = 0) = sin^ [(pd\ 

(17) 

The detection of the spin polarization is affected by the photon polarization. We have the 
relation: 

PL2(^ - <Pd ,) cos2[0] = PLi(^, (f^d) (18) 

In hgure 5 shows the polarizations for A = Cc = 0. Pl^i{6d = f,0d) represents the p 
polarization and P^= 2 (^rf = f)0rf) is the result for the p polarization. 
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p' A=0,e;.=0.] 



FIG. 5: P|(0rf = 4)d) for A = ec = 0. The thick line represents the p polarization and the 

thin line represents the s polarization 

[^=0,^^,<9j =;r/4,e=0, A=0,€c=0.05] 



FIG. 6: Ps=i(0(i = f, 4>d) for Cc = 0.05 at the fixed energy. 

For a finite gap A the polarization is suppressed by the energy factor (sin[/?[e, A)])?. 

For warping figure 6 shows the warping caused by the angle 0^, (sin[/d[e, ec)]) 2 . 

V-Conclusions 

To conclude, a model for computing the photoelectrons intensity and polarization from 
the surface of a topological insulator based on Green’s functions has been introduced. We 
show that the polarization of photoelectrons depends on the laser light polarization in qual¬ 
itative agreement with experimental results [5] . The photoelectrons polarization is modified 
by the incoming photon polarization. This results hold also in the presence of a Zeeman 
gap or warping. For the Zeeman gap the polarization and intensity is suppressed closed to 
the Fermi energy. For warping the polarization and the intensity oscillates with the warping 
angle (pd allowing to identify the spin texture. 

The calculation is based on the tunneling amplitude of the surface electrons into the 
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vacuum .This amplitude can be estimated from the inverted and vacuum TI gap. We com¬ 
pute the number of the emitted photoelectrons. This number depends tunneling amplitude, 
chemical potential,incoming number of photons and weakly dependent on the location of the 
detector. Our calculation indicate that for r = 0.1, /i = O.leV, KpL = 10^ and W = 10^^ 
we obtain that the maximum number of photoelectrons is approximately 2. 

The calculation ignores the bulk electrons, therefore the intensity computed can not be 
compared with the experimental results for energy e — /x < —fi (the bottom of the surface 
states are at Cq = fiv\K\ = 0 , in figures 1 — 4 this corresponds to e — /i = —O.leV ) 
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